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Abstract. In this note we summarize the results from a longer article on obtaining the QCD 
phase diagram as a function of the temperature and chemical potential at large Nc and large 
^ |. Nf in the weak coupling limit A — >■ 0, and the strong coupling limit A oo. The weak 

O ' coupling phase diagram is obtained from the Polyakov line order parameter, and the quark 

^ ^ number, calculated using 1-loop perturbation theory for QCD formulated on 5^x5*^. The strong 

coupling phase diagram is obtained from the same observables calculated at leading order in the 
lattice strong coupling and hopping parameter expansions. We show that the matrix models 
^ , in these two limits agree at temperatures and chemical potentials which are not too high, such 

' that observables in the strongly-coupled theory can be obtained from the observables in the 

' weakly-coupled theory, and vice versa, using a simple transformation of variables. 
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en ' 1. Introduction 

QCD at non-zero chemical potential provides a description of systems at large densities, yet it 
is a description which is currently not directly accessible since it occurs at strong coupling and 
^ ' the non-zero chemical potential leads to a complex action, giving rise to the well-known sign 

^ ■ problem. What this means is that the conventional techniques of studying finite temperature 

QCD: conventional lattice simulations and ordinary perturbation theory, are not applicable. In 
this note we consider two nonconventional perturbative techniques which allow us to calculate 
the partition function and related observables such as the Polyakov line and quark number, and 
from these to map out the phase diagram in an otherwise inaccessible range of temperatures 
and chemical potentials. This is an executive summary of our longer paper pj. 

The two perturbative techniques we employ are 1) weakly-coupled QCD from continuum 1- 
loop perturbation theory on x S^, and 2) strongly coupled lattice QCD with heavy quarks at 
leading order in a strong coupling, and hopping parameter expansion. In both cases we perform 
the calculations in the Veneziano large Nc limit, where Nf oo, while remains fixed. 
This leads to advantages in both cases (for recent reviews on the progress towards understanding 
gauge theories and large Nc see [1, 01 ) • In the weakly coupled theory on x S^, the large Nc 
limit is required to have sharp phase transitions, since the calculation is only valid in very small 
volumes, such that R <^ Ag^^, where R is the radius of the S^. For the lattice strong coupling 
expansion, the spatial volume is large, but large Nc factorization and translational invariance 



lead to a simplification of the action in that terms which include correlations between different 
lattice sites drop out. In both cases the large Nc limit allows for a description of the theory 
in terms of the distribution of the Polyakov line eigenvalues such that the theory reduces to 
an analytically solvable matrix model. For temperatures and chemical potentials which are 
not too high we will show that there is an exact correspondence of the matrix models of the 
weakly-coupled and strongly-coupled theory, under a simple change of parameters. 

2. QCD on x with A — ^ vs. lattice QCD with heavy quarks as A — > oo 

The action for continuum 1-loop QCD with constant was derived in for theories with 
a matter content of scalars, vectors, and/or fermions. For QCD with Nf quarks of mass m, 
chemical potential fj,, and at temperature T = ^, the action in terms of the Polyakov line 

observable /)„ = ^ Tr ^e" io" <^t^o{x) ^ ^ j^fj^ e™^» takes the form 0, 

oo 



Ss^xS^ — Svdm = — —ZvnPnP-n 
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n=l 



(1) 



n 

n=l 



where Sydm is the contribution from the Vandermonde determinant, and Zt,„, Zfn refer to the 
single particle partition functions for vectors and fermions. 
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^^„ = 2j;/(/ + l)e-"^V('+^)'+-^^'. (3) 
1=1 

The action for large Nc-, large Nf lattice QCD in terms of the Polyakov line Vl^(x) = 
Tr P3^Q~^ C/fj, at leading order in the strong coupling and hopping parameter expansion is 
given by [1,13] 

Slat - Svdm = -JDY, [{W)W\x) + {W^)W{x) - {W){W^) 

r 1 (4) 

X 

[ &!A\ fnr inverse cminlincr Ri.-i = 



where J = 2 {^-^j^j for inverse coupling (3iat = and number of temporal slices A',-, and 

h = 2^K^^ is the hopping parameter with k = a^^i_^£) for lattice spacing a and number of 
spatial dimensions D. 

The actions in ([T]) and (HD appear fairly similar with the exception of the sum over n in 
the former, and the sum over x in the latter (the term {W)W^{x) + {W^)W{x) — {W){W^) in 
dH compared to piP-i in ([T]) leads to the same equations of motion). The lack of terms with 
correlations between different lattice sites in Q means that observables of the form {F{W, W'^)) 
will undergo large cancellations such that 

{F{W,W^)) = / J|dVF(x)e-^[^(^)'^^(^)] ^F[VF(x'),VF"^(x')] , 

Co) 

_ /dTye-^(^'^^)F(H^,T^t) ^ ' 

/dVFe-^(W^'^^) 



Therefore, when it is possible to truncate the sum over n in ([T]) to the n = 1 term there is an 
exact correspondence of matrix models resulting from ([1]) and (jH) under the transformations 



Truncation to the 
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z^i ^ JD 



(6) 



1 contribution in ([T]) is valid when the temperature is not too 



and the chemical potential is not too high {fi eji 



1 rnR-^Q^' "^^^^ region of validity includes the line of transitions extending 
from the temperature-axis to the chemical potential-axis, but corrections from Zfn for n > 1 
would need to be included to go to higher chemical potentials, as done in [5|]. 
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3. Large N^, formalism 

To calculate observables in the large limit we adopt the methods in which adapts 
the Gross- Witten-Wadia [1, 0] formalism to handle theories with complex actions. Consider 



a contour C with the Polyakov line eigenvalues Zj 
defined according to the map 



e i distributed along it with density q{z) 
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(7) 



where = vr in the confined phase, such that the contour is closed, and ijj < tt in the deconfined 
phase, such that the contour is open and the distribution of eigenvalues z(s) has a gap. This is 
illustrated in Figure [T] as the deconfinement transition is crossed for /i 7^ 0. 

The blue curve is a distribu- 
tion z{s), obtained from i J ds = 
f dzQ{z), where g{z) is obtained 
in the confined phase at a point 
close to the deconfinement tran- 
sition. Notice that the eigenvalue 
distribution lies significantly away 
from the unit circle (green dotted 
curve) in the —z direction. The 
red curve is the distribution z{s) 
obtained from q{z) in the decon- 
fined (gapped) phase, at a point 
just past the deconfinement tran- 
sition. These distributions corre- 
spond to configurations with com- 
plexified gauge fields (the 6j are 
complex), which is necessary to 
obtain the correct stationary so- 
lution since the action is complex. 

Another requirement for ob- 
taining the correct saddle point solutions is applying the SU{Nc) constraint. Without this 




Figure 1. Distribution z{s) of the Polyakov line eigenvalues 
in the confined phase (blue, fiR = 0.74), and the deconfined 
phase (red, ^lR = 0.76), for TR = 0.3, ^ = 1, mi? = 0. 



constraint taking ^ 7^ trivially shifts ^0 by a constant such that the free energy is independent 
of and the quark number is always zero. The SU{Nc) constraint is incorporated by adding an 
appropriate term to the action with a Lagrange multiplier M, 

S + iAfNc 9i , (8) 

i=l 

where Af = -wiNq is the effective quark number in the large Nc limit [5]. The density q{z) is 

obtained by solving the equation of motion from ^ = 0, which, using ([8]) with ([1]) or takes 
the form 

^ / =^.Q{z')^ = a-^z-a^z-^-M, (9) 

J ZTT^ Z — Z 

where *P indicates principal value and a±\ = Zyip±i + ^Zfie^^^. The procedure for solving 
the equation of motion depends on whether the contour C is open, as in the deconfined phase, 
or closed, as in the confined phase. 



3.1. Confined (ungapped) phase 

When C is considered to be a closed contour, as in the confined phase, the equation of motion ^ 
can be solved for g{z) using Cauchy's theorem. In the confined phase it is sufficient to consider 
the Fourier expansion of the density 

00 

q{z)= PnZ""'^ (10) 

n=— 00 

and solve for the pn- The density is constrained to satisfy the identity constraint 

tfT > / —g{z) = i, (11) 

as well as the SU (Nc) constraint 

r d7 

Y0. = > / -^g{z) log(z) = . (12) 

The free energy and other relevant observables can then be calculated by plugging in the 
stationary point solutions obtained for the p„. 



3.2. Deconfined (gapped) phase 

When C is considered to be an open arc, as in the deconfined phase, then the equation of motion 
dHl must be solved by defining a resolvent and solving the Plemelj formulae. The resolvent is 
defined from the singular integral contribution in the equation of motion 



Following [3, ISl] the contour C along which the eigenvalues are distributed is defined as a square 
root branch cut. The resolvent can then be evaluated using singular integral techniques for open 
arc contours to obtain [1|] 



(a-iz^-Qiz ^) + V z"^ + - 2rxz [air ^z ^ + a_i) , (14) 



where the endpoints of the arc occur at radius r and angle with x = cosip. Observables 
can then be calculated using 

where T is defined as a contour around C, which can then be peeled off to surround the residues 
outside such that Cauchy's theorem can be used. 

4. Results 

Using the methods of the previous sections we calculated the Polyakov lines in the confined 
and deconfined regions and mapped out the phase diagrams [1]. Figure [2] shows the phase 
diagram for the weakly-coupled theory and Figure [3] for the strongly-coupled theory. The phase 
boundaries are determined by comparing the free energies of the gapped and ungapped eigenvalue 
distributions in the regions where both are possible. In the ungapped region the effective quark 
number J\f = and where the gapped distribution is favored M ^ 0. In both phase diagrams 
the order of the transition is at least fifth order at /u = 0, and at least 2nd order for 7^ 0. 
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Figure 2. Phase diagram on x S"^ for 
mi? = 0, ^ = 1. 
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Figure 3. Phase diagram from the strong 
coupling expansion for h = 0.01. 
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